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Abstract 

Various technical developments enlarged the potential of angle-resolved photo emission (ARPES) tremendously dur- 
ing the last one or two decades. In particular improved momentum and energy resolution as well as the use of photon 
energies from few eV up to several keV makes ARPES a rather unique tool to investigate the electronic properties of 
solids and surfaces. Obviously, this rises the need for a corresponding theoretical formalism that allows to accompany 
experimental ARPES studies in an adequate way. As will be demonstrated by several examples this goal could be 
achieved by various recent developments on the basis of the one-step model of photo emission: The spin-orbit in- 
duced Rashba-splitting of Shockley-type surface states is discussed using a fully relativistic description. The impact 
of chemical disorder within surface layers can be handled by means of the Coherent Potential Approximation (CPA) 
alloy theory. Calculating phonon properties together with the corresponding electron-phonon self-energy allows a 
direct comparison with features in the ARPES spectra caused by electron-phonon interaction. The same holds for 
the influence of electronic correlation effects. These are accounted for by means of the dynamical mean field theory 
(DMFT) that removes the most serious short comings of standard calculations based on the standard local density 
approximation (LDA). The combination of this approach with the CPA allows the investigation of correlated transi- 
tion metal alloys. Finally, accounting for the photon momentum and going beyond the single scatter approximation 
for the final state allows to deal quantitatively with ARPES in the high-energy regime (HAXPES) that reduces the 
influence of the surface on the spectra and probing primarily the bulk electronic structure this way. Corresponding 
calculations of ARPES spectra, however, have to deal with thermal vibrations in an adequate way. For this, a new 
scheme is suggested that makes use of the CPA. 
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1. Introduction 

The spectrum of one-particle excitations of a sys- 
tem of correlated electrons in a solid is a fundamental 
question in condensed-matter physics. The theoretical 
understanding of the excitation spectrum poses a long- 
standing and not yet generally solved problem. Within 
the independent-electron approximation the spectrum is 
simply given in terms of the one-particle eigenenergies 
of the Hamiltonian. Analogously, it is widely accepted 
to interpret a measured photo emission spectrum by re- 
ferring to the results of band-structure calculations that 
are based on density functional theory (DFT) and the lo- 
cal density approximation (LSDA) 01 [US]. Such an in- 
terpretation is questionable since there is actually no di- 
rect correspondence between the Kohn-Sham eigenen- 
ergies and the one-particle excitations of the system 



Jl, 01 • For an in principle correct description of the ex- 
citation energies, the local LSDA exchange-correlation 
potential has to be supplemented by the non-local, com- 
plex and energy-dependent self-energy which leads to 
the Dyson equation J3I IH instead of the Schrodinger- 
or Dirac-type equation in the non- or fully relativistic 
Kohn-Sham scheme. 

To overcome these deficiencies we have developed a 
generalized approach by accounting properly for elec- 
tronic correlations beyond the LSDA. For this pur- 
pose a general non-local, site-diagonal, complex and 
energy-dependent self-energy Y. DMFT |0, [jj, has 
been included self-consistently in the fully relativis- 
tic Korringa-Kohn-Rostoker multiple scattering theory 
(SPRKKR) [9]. Within this scheme the self-energy is 
determined by a self-consistent dynamical mean field 
theory (DMFT) calculation. This method is straight- 
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forwardly applicable for electronic structure calcula- 
tions of 3D- and 2D-systems with perfect lateral transla- 
tional invariance and arbitrary number of atoms per unit 
cell. Furthermore, the calculational scheme has been 
extended to the case of disordered alloys |9|], to include 
many-body correlation effects in the electronic structure 
and photo emission calculations of these class of com- 
plex compound materials. 

Provided that the self-energy is known, one can de- 
duce a PES/IPE "raw spectrum" from the solution of 
the Dyson equation. To achieve a reliable interpreta- 
tion of experiments, however, it is inevitable to deal 
with so-called matrix-element effects which consider- 
ably modify and distort the raw spectrum. Above all, the 
wave-vector and energy dependence of the transition- 
matrix elements has to be accounted for. These depen- 
dencies are known to be important and actually can- 
not be neglected. They result from strong multiple- 
scattering processes which dominate the electron dy- 
namics in the low-energy regime of typically 1-200 eV 
ifioll . The transition-matrix elements also include the 
effects of selection rules which are not accounted for in 
the raw spectrum. Strictly speaking, it can be stated that 
the main task of a theory of photo emission is to close 
the gap between the raw spectrum obtained by LSDA 
or LSDA+DMFT electronic-structure calculations and 
the experiment. The most successful theoretical ap- 
proach is the so-called one-step model of photo emis- 
sion as originally proposed by Pendry and co-workers 
iTToL [0 , 12 ] . In the following a short overwiew will 
be given on the recent extensions of the one-step model 
which are connected with correlation effects, disordered 
alloys, phonon related features and high excitation en- 
ergies far beyond the ultraviolet regime. 

2. Recent developments in the one step model of 
photo emission 

The main idea of the one-step model is to describe 
the actual excitation process, the transport of the photo- 
electron to the crystal surface as well as the escape into 
the vacuum 111 311 as a single quantum-mechanically co- 
herent process including all multiple-scattering events. 
Within this model self-energy corrections, which give 
rise to damping in the quasi-particle spectrum, are prop- 
erly included in both the initial and the final states. This 
for example allows for transitions into evanescent band 
gap states decaying exponentially into the solid. Sim- 
ilarly the assumption of a finite lifetime for the initial 
states gives the opportunity to calculate photo emission 
intensities from surface states and resonances. Treat- 
ing the initial and final states within the fully relativistic 



version of layer Korringa-Kohn-Rostoker (KKR) theory 
01411 . it is a straight forward task to describe complex 
layered structures like thin films and multilayers within 
the photo emission theory. Furthermore, the surface de- 
scribed by a barrier potential can be easily included into 
the multiple-scattering formalism as an additional layer. 
A realistic surface barrier model which shows the cor- 
rect asymptotic behavior has been introduced, for exam- 
ple, by Rundgren and Malmstrom 1 15]. 

We start our considerations by a discussion of 
Pendry 's formula for the photocurrent which defines the 
one-step model of PES 111 111 : 



oc Im (E f , k|||GjAG|A t G2 \E f , ky> 



(1) 



The expression can be derived from Fermi's golden rule 
for the transition probability per unit time |4], with 7 PES 
denoting the elastic part of the photocurrent. Vertex 
renormalizations are neglected. This excludes inelastic 
energy losses and corresponding quantum-mechanical 
interference terms 1611 . Furthermore, the inter- 

action of the outgoing photoelectron with the rest sys- 
tem is not taken into account. This "sudden approx- 
imation" is expected to be justified for not too small 
photon energies. We consider an energy-, angle- and 
spin-resolved photo emission experiment. The state of 
the photoelectron at the detector is written as \Ef, ky), 
where ky is the component of the wave vector parallel 
to the surface, and E y is the kinetic energy of the pho- 
toelectron. The spin character of the photoelectron is 
implicitly included in |£/,k||) which is understood as 
a four-component Dirac spinor. The advanced Green 
function in Eq. ([TJ characterizes the scattering prop- 
erties of the material at the final-state energy E% = Ef. 
Via ^f) - G^\Ef,k\\) all multiple-scattering correc- 
tions are formally included. For an appropriate descrip- 
tion of the photo emission process we have to model the 
correct asymptotic behavior of *F /(r) beyond the crystal 
surface, introducing a single outgoing plane wave char- 
acterized by Ef and ky . Furthermore, the damping of the 
final state due to the imaginary part of the inner potential 
iVoiiEi) must be taken into account. We thus construct 
the final state within spin-polarized low-energy electron 
diffraction (SPLEED) theory considering a single plane 
wave \Ef, ky) advancing onto the crystal surface. Using 
the standard layer- KKR method 111411 generalized for the 
relativistic case 1 17, 3, we first obtain the SPLEED 
state -r*P/(r). The final state is then given as the time- 
reversed SPLEED state (T = -icr y K is the relativis- 
tic time inversion operator). Many-body effects are in- 
cluded phenomenologically in the SPLEED calculation, 
by using a parametrized, weakly energy-dependent and 
complex inner potential Vo(£2) = Vo r (E2) + iVoi(E2) 
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ifioll . This generalized inner potential accounts for in- 

I — I 

elastic corrections to the elastic photocurrent [4] as well 
as the actual (real) inner potential, which serves as a 
reference energy inside the solid with respect to the 
vacuum level II 1 911 . Due to the finite imaginary part 
Vo;(£ , 2), the flux of elastically scattered electrons is con- 
tinuously reduced, and thus the amplitude of the high- 
energy wave field ^/(v) can be neglected beyond a cer- 
tain distance from the surface. 

The practical calculation starts with the Dirac Hamil- 
tonian /zlsda which one has to consider in the frame- 
work of DFT 1I20L 21]. In atomic units (h — m — e — 



with 



1 , c — 1 37.036) one has: 

/jLSDA(r) = -icaV +/3c 2 -c 2 + V hSDA (r) + /3(rB hSDA (r) . 

(2) 

VlsdaW denotes the (effective) spin-independent po- 
tential, and BlsdaW is the (effective) magnetic field 
given by B22I1 : 



VLSDA(r) = ^ SDA W + Vi SDA (r)), 



(3) 
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l(Vh nA (r) - V^{r))e B . (4) 



2 v LSDA V 

The constant unit vector e# determines the spatial direc- 
tion of the (uniform) magnetization as well as the spin 
quantization axis. /3 denotes the usual 4x4 Dirac matrix 
with the nonzero diagonal elements f}\\ = P22 = 1 and 
/?33 = Paa - — 1, and the vector a is given by its com- 
ponents at = <j x ® eric (k = x, y, z) in terms of the 2x2 
Pauli-matrices at- 

The "low-energy" propagator G\ in Eq. (Q~|), i. e. the 
one-electron retarded Green function for the initial state 
in the operator representation, yields the "raw spec- 
trum". It is directly related to the "bare" photocur- 
rent and thereby represents the central physical quan- 
tity within the one-step model. G\ = G\(Ei) is to be 
evaluated at the initial-state energy E\ = Ef — co, where 
u> is the photon energy. In the relativistic case G\ is 
described by a 4 x 4 Green matrix which has to be ob- 
tained for a semi-infinite stack of layers. This quantity 
is defined by the following equation: 



[Ei - /i LS DA(r)]Gf(r, r',£,) =6(r-r') 



(5) 



In order to account for strong electronic correlations 
beyond the LSDA-scheme one has to introduce a non- 
local, energy and spin-dependent potential V. This quan- 
tity can be defined in the following way: 

V(r, r', E) = tf(r - r') (V LSD A(r) + /3crB LSDA (r)) 
+ Z w (r, r', E) + p&l^Xr, r', E) , (6) 



S (y) (r, r', E) = ^0(r, r', E) + 2 j (r, r', E)) (7) 

and 

S (B) (r, r', E) = ^(£ T (r, r', E) - I x (r, r', E)) e B . (8) 

The resulting integro-differential equation for the 
initial-state one-electron retarded Green function takes 
the form: 

\e + icaV -f3c 2 + c 2 )] G|(r, r', E) 
+j V(r, r", E) G[(r", r', E)dr" = S(r - r') . (9) 

According to the LSDA+DMFT approach realized in 
the framework of the fully relativistic Korringa-Kohn- 
Rostoker multiple scattering theory (SPR-KKR) |9J we 
use a self-energy Z® MFT {E) calculated self-consistently 
using dynamical mean-field theory (DMFT) iSSHtl- 
The explicit form in relativistic notation is given by: 

V AA ,(r,r',£) = (V(r) + &B(r)) 6(r - r') 

+ lP h f ir T (E)5a5 ll2 . (10) 

The spin-orbit quantum number k and magnetic quan- 
tum number /u were combined in the symbol A = (K,fj). 

For disordered alloys the coherent potential approx- 
imation (CPA) scheme can straightforwardly be ex- 
tended to include the many-body correlation effects for 
disordered alloys [9]. According to the LSDA+DMFT 
approach realized in the framework of the fully relativis- 
tic SPR-KKR multiple scattering theory we use a self- 
energy T.® MFT (E) calculated self-consistently using dy- 
namical mean field theory la, |7I, la] f° r trie component a 
of an alloy. Within the KKR approach the local multi- 
orbital and energy dependent self-energy E^ MFT (E) is 
directly included into the single-site matrix t a for the 
component a, respectively by solving the corresponding 
Dirac Eq. ([9j- Consequently, all the relevant physical 
quantities connected with the Green's function and used 
for photo emission calculations, as for example, the ma- 
trix elements or scattering matrices like t CPA contain the 
electronic correlations beyond the LSDA scheme. A de- 
tailed description of the generalizated one-step model 
for disordered magnetic alloys can be found in llzil . 
Here we present the expression for the CPA photocur- 
rent only. Following Durham et al. 1124 12511 we have 
to perform the configurational average for the various 
contributions to the photocurrent: 

</ PES (e / ,k||)> = </ a ( e/ ,k||)> + </ m (e / ,k||)> 
+ </ s (e / ,k||)> + </ mc (e / ,k||)> 

(11) 
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When dealing with disorder in the alloys, an additional 
contribution 7 lnc , the so called "incoherent" term ap- 
pears. This contribution to the alloy photocurrent ap- 
pears because the spectral function of an disordered al- 
loy f26ll is a non single-site quantity. In fact this con- 
tribution is closely connected with the presence of the 
irregular wave functions well-known from the spherical 
representation of the Green function G|. 

Within the next section we present some selected re- 
sults from photo emission investigations on various sys- 
tems to demonstrate the great flexibility of the method 
described above. We start with a photo emission study 
on MgO/Fe(001) performed for photon energies in the 
hard X-ray regime. 

3. High energy photo emission spectroscopy 



model has been extended to 
Starting with a new 



28] 



Recently, the one-ste 
the soft X-ray regime 
ansatz for the final state which generalizes the so called 
single-scatter approach, it becomes possible to account 
for the lattice for the complete photo emission calcula- 
tion. This method is applicable also for higher photon 
energies ranging from 1 keV up to more or less than 
10 keV. Increasing the photon energy means increasing 
the bulk sensitivity of the corresponding photo emission 
data. We demonstrate this effect by comparing spectra 
obtained by photoelectron excitation with 1 keV and 6 
keV radiation from the clean Fe(001) surface and from 
the overlayer system MgO/Fe(100) with 8 ML MgO on 
Fe. Fig. la presents the spectrum for the clean Fe(001) 
surface and a photon energy of 1 keV. As expected for 
this photon energy regime, the bulk sensitivity is en- 
hanced and the surface emission is reduced to a negligi- 
ble extent. Also the relative intensity fraction of the sp- 
bands is obviously increased. Going to a much higher 
photon energy of 6 keV the d-band intensity is strongly 
reduced when compared with sp-band related features. 
This is shown in Fig. lb. In a next step we put 8 ML of 
MgO on the Fe(001) surface and repeat our calculations 
for both photon energies. This is shown in the lower 
panel of Fig. 1. At 1 keV mostly the MgO bulk band- 
structure is visible in Fig. lc. This is due to the thickness 
of the MgO film which consists of 8 ML. Increasing 
the photon energy to 6 keV we expect due to the much 
larger mean free path length of the photo electron that 
Fe-related features will recover, namely the Fe sp-states. 
This is clearly demonstrated by Fig. Id which represents 
the corresponding photo emission spectrum calculated 
for the 8ML MgO/Fe(001) system at a photon energy of 
6 keV. Besides the effect that more than one Brillouin 
zone is visible for a fixed escape angle regime due to 




Figure 1 : (Color online) Photoemission spectra calculated for clean 
Fe(001) (upper panel, Fig. la and lb). Left side shows the intensity 
distribution obtained for 1 keV, right side represents the correspond- 
ing spectrum at 6 keV. Lower panel shows theoretical spectra for the 
overlayer system 8MgO/Fe(001) at 1 keV (left side, Fig. lc) and 6 
keV (right side, Fig. Id). 



the very high photon energy it is undoubtly observable 
that nearly all MgO related features are vanished in the 
intensity plot. This, for example, should give access to 
buried interfaces. 



4. Treatment of correlation effects within dynamical 
mean field theory (DMFT) 

The following examples concerns the ferromagnetic 
transition metal systems Ni and Fe as prototype mate- 
rials to study electronic correlations and magnetism be- 
yond the LSDA scheme. In Fig. 2 we present a com- 
parison between experimental photo emission data i29ll 
and calculated spectra using different theoretical ap- 
proaches 113011 . In the upper row spin-integrated ARUPS 
measurements from Ni(Oll) along FY for different an- 
gles of emission are shown. The dotted lines represent 
the experimental data, whereas the solid lines denote 
the single-particle approach to the measured spectral 
function. Obviously, the LSDA-based calculation com- 
pletely fails to describe the experimental data. The ener- 
getic positions of the theoretical peaks deviate strongly 
from the measured ones. Furthermore, the compli- 
cated intensity distributions that appear for higher an- 
gles of emission could not be reproduced by the LSDA- 
based calculations at all. In contrast, the non self- 
consistent quasi-particle calculation provides a signif- 
icant improvement when compared to the measured 
spectra. For the complete variety of emission angles 
the energetic peak positions coincidence with the ex- 
periment within about 0. 1 e V. Only the overall shape of 
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Figure 2: (color online) Spin-integrated ARUPS spectra from Ni(01 1) 
along T7 for three different angles of emission. Upper row: com- 
parison between LSDA-based calculation and experiment [29(1; mid- 
dle row: comparison between experiment and non self-consistent 
quasi-particle calculations neglecting matrix element and surface ef- 
fects 12911 : lower row: spin-integrated LSDA+DMFT spectra includ- 
ing photo emission matrix elements (this work). Theory: solid red 
line, experiment: black dots. 
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Figure 3: (a) Experimental spin-integrated photo emission spectra of 
the Fe(110) surface measured with p-polarization in normal emission 
along the TN direction of the bulk Brillouin zone. The curves are la- 
beled by the wave vectors in units of TN=1.55 A~'. (b) Correspond- 
ing one-step model calculations based on the LSDA+DMFT method 
which include correlations, matrix elements and surface effects. 



the measured spectral intensities deviate from the calcu- 
lations because of the neglection of multiple scattering 
and surface-related effects. In the experiment the differ- 
ent peaks seem to be more broadened and the spectral 
weight especially for nearly normal emission is shifted 
by about 0.1 eV to higher binding energies. In addi- 
tion it seems that for very high emission angles like 60° 
an even more complicated peak structure is hidden due 
to limited experimental resolution. An additional spin- 
analysis is therefore highly desirable for these experi- 
ments. 

Within our work we could go far beyond previous 
theoretical studies by combining a sophisticated many- 
body approach as the self-consistent LSDA+DMFT 
method with the one-step based calculation of the corre- 
sponding spectral function. The self-energy within the 
DMFT has been applied only for d-states and can be 
calculated in terms of two parameters - the averaged 
screened Coulomb interaction U and the exchange in- 
teraction J. The screening of the exchange interaction 
is usually small and the value of J can be calculated 
directly and is approximately equal to 0.9 eV for all 3d 
elements. This value has been adopted for all of our cal- 
culations presented here. For U we used a value of 3 eV. 
However, our test calculations show that a somewhat 



different choice of U does not substantially changes 
features and trends in the calculated spectra. The in- 
tensity distributions resulting from the corresponding 
photo emission calculation are shown in the lower row 
of Fig. 2. A first inspection reveals very satisfying quan- 
titative agreement between experiment and theory for 
all emission angles. Let us concentrate first on the ex- 
citation spectrum calculated for = 5°. The spin- 
integrated spectrum exhibits a pronounced double-peak 
structure with binding energies of 0.1 eV and 0.3 eV. 
The second peak is slightly reduced in intensity which is 
also in accordance with the experimental findings. Fur- 
thermore, the width of the spectral distribution is quanti- 
tatively reproduced. The calculated binding energies are 
related to the real part of the self-energy that corrects the 
peak positions due to dynamical renormalization proce- 
dure of the quasiparticles which is missing in a typical 
LSDA-based calculations. The relative intensities and 
the widths of the different peaks, on the other hand, must 
be attributed to the matrix-element effects which enter 
our calculations from the very beginning via the one- 
step model of photo emission. The double-peak struc- 
ture originates from excitation of the spin-split d-bands 
in combination with a significant amount of surface- 
state emission ll3lll . The two spectra calculated for high 
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angles of emission show the more broadened spectral 
distributions observable from the experimental data. An 
explanation can be given in terms of matrix-element ef- 
fects, due to the dominating dipole selection rules. The 
spin-resolved spectra reveal a variety of d-band excita- 
tions in both spin channels, which in consequence lead 
to the complicated shape of the spectral distributions 
hardly to be identified in the spin-integrated mode. 

Our second example within this section concerns a 
spectroscopic study on ferromagnetic Fe 13211 . Figures 
3a and 3b display a comparison between spin-integrated 
ARPES data and theoretical LSDA+DMFT based one- 
step photo emission calculations of Fe(llO) along the 
TN direction of the bulk Brillouin zone (BZ) with p- 
polarized radiation. In our LSDA+DMFT investigation 
underlying the ARPES calculations we use for the av- 
eraged on-site Coulomb interaction U a value 17=1.5 
eV which lies between the experimental value U *1 eV 
13311 and a value U «2 eV derived from theoretical stud- 



ies II34L 13511 . The k values were calculated from the used 
photon energies ranging from 25 to 100 eV. Near the F 
point (k~0.06 FN), the intense peak close to the Fermi 
level corresponds to a sj 3 minority surface resonance, 

as indicated on top of Fig. 3. Experimentally, its £3 bulk 
component crosses the Fermi level at k ~0.33 FN, lead- 
ing to a reversal of the measured spin-polarization and 
to a strong reduction of the intensity at k =0.68 FN in 
the minority channel. The peak at the binding energy 
BE~0.7 eV, visible mainly for p-polarization in a large 
range of wave vectors between F and N, can be assigned 
to almost degenerate Zj . bulk-like majority states. A 
feature at BE~1.1 eV dominates the spectrum at the T- 
point. Depending on the polarization its degenerate Ej 
states form a shoulder around the same BE. The broad 
feature around 2.2 eV, visible at various Appoints, but not 
at the N-point, is related to a majority E] 3 surface state. 
Around the N-point (0.76< k <1.0) and at BE>3 eV we 
observe a Sj band having strong sp character. The pro- 
nounced difference between its theoretical and exper- 
imental intensity distributions can be attributed to the 
fact that in the present calculations only local Coulomb 
repulsion between d electrons is considered, without ad- 
ditional lifetime effects for the sp bands. Finally, we 
notice that the background intensity of the spectrum at 
k-0.66 FN, corresponding to a photon energy of 55 eV, 
is strongly increasing by the appearance of the Fe 3p 
resonance. 

In conclusion, we have presented spectral function 
calculations of ferromagnetic Ni and Fe, which coher- 
ently combine an improved description of electronic 
correlations with multiple-scattering, surface emission, 
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Figure 4: ARUPS spectra taken from the Ni v Pdi_ v (00I) alloy sur- 
faces as a function of the concentration x for a fixed photon energy of 
/iv=40.0 eV along TX in normal emission. Experimental data shown 
in the left panel calculated spectra presented in the right panel. De- 
pending on the concentration x a pronounced shift in spectral weight 
towards the Fermi level is visible. 



dipole selection rules and other matrix-element related 
effects that lead to a modification of the relative photo 
emission intensities. As has been demonstrated, this 
approach allows on the one hand side a detailed and 
reliable interpretation of high-resolution angle-resolved 
photo emission spectra of 3d-ferromagnets. On the 
other hand, it also allows for a very stringent test of new 
developments in the field of DMFT and similar many- 
body techniques. 

5. Treatment of disordered alloys via the coherent 
potential approximation (CPA) 

In this section we want to study the alloying effect in 
combination with electronic correlations. Fig. 4 shows a 
series of spectra as a function of the concentration x cal- 
culated for a photon energy hv-40 e V with linear polar- 
ized light. The experimental data have been presented 
in the left panel and the corresponding LSDA+DMFT- 
based photo emission calculations are presented in the 
right one. Our theoretical analysis shows that starting 
from the pure Ni case, the agreement is fully quanti- 
tative with deviations less than 0. 1 e V binding energy, 
as expected. Going to the Nin .soPdo.20 alloy the agree- 
ment is on the same level of accuracy concerning the 
range of binding energies between the Fermi energy and 
2 eV binding energy. Inspecting the density of states 
(DOS) for the Nio.80Pdo.20 alloy this fact becomes ex- 
plainable, because this energy interval represents the 
Ni-dominated region. The Pd-states start to appear at 
about 2 eV below Ep besides of the small dip at the 
Fermi level. For higher binding energies the agreement 
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Figure 5: ARUPS spectra taken from the Ni, Aui_ v (l 1 1) alloy sur- 
faces as a function of the concentration x for a fixed photon energy of 
hv=2\.0 eV along YX. The upper panel shows the majority surface 
state, the lower panel shows the minority surface state with increasing 
magnetic moment from left to right. 



is also very good, although a bit more structure is ob- 
servable in the theory, especially around 3.5 eV. An ex- 
planation for this behavior can be found in terms of life- 
time effects, but it should be mentioned here that the 
experimental intensity background consisting of sec- 
ondary electrons was not considered in the theoretical 
analysis. From the Nio.70Pdo.30 alloy system it becomes 
clearly visible that the deviation between theory and ex- 
periment is mainly introduced for increasing concentra- 
tion of Pd. The same argument holds for the Nio.50Pdo.50 
and Nio.30Pdo.70 alloys shown next in the series. In ad- 
dition, the spectra of Nio.30Pdo.70 reveal some deviations 
near the Fermi level. Also, the spectral intensity of the 
Ni surface resonance, that appears at about 0.5 eV bind- 
ing energy is underestimated in the calculation when 
compared to the experiment. 

Our spectroscopical analysis has clearly demon- 
strated that the electronic properties of the NiJPdi_ x 
alloy system depend very sensitively on the interplay 
of alloying and electronic correlation. A description 
within the LSDA approach in combination with a CPA 
method results in a qualitative description of the elec- 
tronic structure of Ni x Pdi_ x only. 

The second example concerns more surface sensitive 
features, namely the spin-orbit-splitting observable for 
surface states from high Z-materials like Au (3j|] ■ Here 
we want to study the interplay between spin-orbit in- 
teraction and correlation using as a prototype material 
the magnetic Au v Nii_ v alloy. We are able to tune the 
magnetic moment of the alloy as a function of the con- 
centration x and in consequence we can study both in- 
teractions on the same level of accuracy. Fig. 5 shows a 



series of spin-resolved photo emission spectra from the 
Au v Nii_ x (l 1 l)-surface. For a very low Ni concentration 
of about 5% the magnetic moment is nearly zero and we 
observe the pure spin-orbit split Shockely surface state, 
well known from literature 13611 . With increasing the Ni 
concentration the magnetic moment increases and the 
magnetic exchange interaction starts to compete with 
the spin-orbit splitting. Exchange causes a splitting in 
energy, whereas spin-orbit interaction causes a splitting 
in k||. For Ni concentrations which lead to an exchange 
splitting comparable with the spin-orbit one we find the 
interesting result that the two parabolas seem to shift 
both in energy and ky with the additional effect of struc- 
tural distortion. For larger magnetic moments the ex- 
change interaction starts to win against the spin-orbit 
interaction and one can observe as usual the exchange- 
split Shockely surface state which has to be expected on 
the Ni(lll) surface. 

These two examples may illustrate that the use of 
the CPA alloy theory self-consistently combined with 
a DMFT+LSDA approach serves as a powerful tool for 
electronic structure calculations, whereas the applica- 
tion of the fully relativistic one-step model of photo 
emission, which takes into account chemical disorder 
and electronic correlation on equal footing guarantees a 
quantitative analysis of the corresponding spectroscopi- 
cal data. 

6. Effects of efectron-phonon interaction in angie- 
resolved photo emission 

Nowadays high resolution photo emission measure- 
ments allow to investigate in very detail the electronic 
structure of materials close to the Fermi level. This 
gives access to energy bands modified by electron- 
phonon interaction, and therefore serves as a very 
important spectroscopical tool to characterize trans- 
port properties or, for ex amp le, prop erties of super- 
conducting materials 11371 l38l l39l 14011 . Modifications 
of the electronic structure generated by electron-phonon 
interaction are typically not accounted for in LSDA- 
based band-structure calculations. The theoretical treat- 
ment of electron-phonon effects requires the use of 
many-body techniques which originally had been devel- 
oped to give a quantitative description of strong elec- 
tronic correlation effects. The self-energy approach of- 
ten applied in many -body calculations can be mapped to 
the case of electron-phonon related correlations. There- 
fore, having the electronic structure obtained by an ab- 
initio calculation, the modifications due to electron- 
phonon interaction can be described via an electron- 
phonon self-energy E(k, E^) I4HI42I1 that appears in the 
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momentum-resolved Green's function operator 
1 



G k (z) = 



z-E k - S(k, E k ) 



(12) 



The electron-phonon self-energy is derived using 
many-body perturbation theory within second order. It 
follows fil: 



ME) 



J 



T. k (E, (jS)a Fk(oS)da> , 



(13) 



where ~Lk{E, cj) represents the self-energy obtained 
within the Einstein model 14 ill . oP'FkicS) denote the 
momentum-resolved Eliashberg function: 

a>F k (a» = J] |^ + " q | 2 x 

q,A,n' 

S(E nM - EpWEn^-n - E F )5{u - Wq) 

(14) 

with the electron-phonon matrix element 

q,A 

= r^— (n ' k K>'> k - q> -( 15 > 



A,n,n' 
^k,k+q 



2M v wt 



Herein, k, n stand for the wave vector and band index, 
respectively, q is the phononic wave vector and A de- 
notes the number of phonon modes. e£' a defines the po- 
larisation of a phonon mode and v numbers the atoms 
within the unit cell. 

To obtain explicitly the electron-phonon self-energy, 
we focus our attention on the Eliashberg function as the 
main quantity determined through all coupled electronic 
and phononic states of the system. This function can be 
reformulated in terms of the Green's function by use of 
the identity (see e.g. Butler (1985)). After some alge- 
bra the Eliashberg function follows as the sum over all 
phonon states: 



a F k (cj) 



^q,A q,A 

y y v ' a vfi Hu-ccfo (i6) 

q,A v,a v q 



with 



T tp,q = f d'rd'r'VVy^G^VVy^G^. (17) 



According to Gyorffy and coworkers B43I1 the Green's 
function G^'* = G^(r, r, Ep ± iO) can be expressed in 
terms of the scattering path operator at the Fermi level 
(t£(2?f + ;0))a'a and is therefore accsessible in terms of 
multiple scattering theory. 







Pb A 
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Figure 6: Upper panel: Calculated Eliashberg function oP'F^O)) for 
Pb(ll())fork = 2(0.17, 1., 0.17). Lower panel: corresponding k- 
dependent self-energy, real part shown on left, imaginary part shown 
on the right. 



The complex self-energy contributes due to electron- 
electron and electron-impurity scattering, as well as due 
to electron-phonon scattering described by the electron- 
phonon self-energy 1* e i-ph- The real and imaginary 
parts of "Lei-ph are responsible for band renormalisations 
around E-p. These quantities are determined by electron- 
phonon scattering induced quasiparticle lifetimes, mea- 
surable in photo emission experiments. To demonstrate 
the effect of electron-phonon interactions on the elec- 
tronic quasi-particle states near Ep we present photo 
emission calculations for Pb, which serves as a proto- 
type system and compare our theoretical results to cor- 
responding experimental data [36]. Pb is of special in- 
terest because of the presence of strong electron-phonon 
interactions which in particular lead to the supercon- 
ducting state at low temperatures I37ll40ll . 

The calculations have been performed for energies 
close to the Fermi level. The wave vectors in the BZ 
(k = ^(0.17, 1 ., 0. 17)) correspond to the experimental 
geometry M4CH1 . Fig.|6]shows a comparison between the 
experimental and theoretical real and imaginary parts of 
the self-energy I. k . Obviously, the agreement is very 
satisfying. 

In general, it is possible to deduce a phononic self- 
energy from measured photo emission spectra. The 
imaginary part of the self-energy can be estimated from 
the full width at half maximum of the energy distribu- 
tion curves fitted by a Lorenzian. The real part is than 
obtained by a Kramers-Kronnig transformation. How- 
ever, in cases where for example several transitions ap- 
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Figure 7: ARUPS spectra taken from the Pb(llO) surface exclud- 
ing (left panel) and including (middle padel) phonons calculated self- 
consistently from first principles. Corresponding experimental data 
are shown on the right panel [ 4411 . 



pear close to the Fermi level or even additional surface 
contributions to the photocurrent take place, this anal- 
ysis might be questionable. Therefore, it seems to be 
quite reasonable to make a direct comparison between 
angle resolved photo emission calculations and corre- 
sponding experimental data. For this purpose we de- 
veloped a scheme that allows to include directly a k- 
dependent self-energy into the one-step model of photo 
emission. In Fig. 7 we present photo emission spec- 
tra calculated for the two cases without (left pannel) 
and with (middle panel) electron-phonon interactions 
and compare them to the correponding measurements 
1 44J] ■ The contour plot has been calculated for Hel light 
along fX More precisely, bands forming the locally 
tubular Fermi surface of Pb(l 10) crosses at about 5° 
off-normal emission along FX. In the middle pannel of 
Fig. 7 a clear renormalisation of the bare band (often 
called kink) at 8meV binding energy is observable. This 
effect appears just where the real part of the self-energy 
reaches its maximum. Futhermore, an additional broad- 
ening of the spectral features due to the imaginary part 
of the phononic self-energy is visible in the calculated 
spectra. This, of course, completes a lifetime analyses 
of spectral features in a quantitative sense because the 
imaginary part of the phononic self-energy can be iden- 
tified as the last contribution that has to be considered 
in addition to impurity scattering and electronic correla- 
tion related broadening mechanisms of photo emission 
intensities. 



7. Conclusions 

We have shown several examples of photo emission 
calculations which reflect the wide range of applicabil- 
ity of the one-step model of photo emission. In its lat- 
est version that is directly combined with the SPRKKR 



program package quantitative spectroscopical investiga- 
tions may range from 10 eV up to 10 keV concerning 
the photon energy. This gives not only direct accsess 
to the UV regime but also to the so called HAXPES- 
regime of angle resolved photo emission. Due to the 
LSDA+DMFT method the range of applicability has 
been greatly enhanced ranging now from simple param- 
agnetic metals to complex layered compounds inhibit- 
ing strong electronic correlations. The recent extent of 
the one-step model to disorderd magnetic alloys round 
up the range of applicability to a variety of complex 
3D- as well as 2D-materials and allows investigations 
on surface related features for low photon energies as 
well as bulk like analyses for higher photon energies. 
As it holds for photocurrent calculations in the HAX- 
PES regime phonon related effects are also very impor- 
tant when analysing the intensity distributions which 
are obtained for energies around the Fermi level. The 
second point has been addressed in the last section in 
which we demonstrated how to introduce in a quati- 
tative sense a k-dependent phononic self-energy that 
allows for spectroscopical calculations from kink-like 
structures. Therefore this ansatz enlarges in combina- 
tion with the LSDA+DMFT method again the range of 
applicability to systems like high T c cuprats or other 
highly correlated classes of materials. 
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